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ABSTRACT 

We propose a differential analog of the notion of integral 
closure of algebraic function fields. We present an algorithm 
for computing the integral closure of the algebra defined 
by a linear differential operator. Our algorithm is a direct 
analog of van Hoeij’s algorithm for computing integral bases 
of algebraic function fields. 
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1. INTRODUCTION 

The notion of integrality is a classical concept in the the¬ 
ory of algebraic field extensions. If R is an integral domain, k 
the quotient field of R, and K = A:(a) an algebraic extension 
of k of degree d, then an element of K is called integral if its 
monic minimal polynomial has coefficients in R. While K 
forms a fc-vector space of dimension d, the set of all integral 
elements of K forms an i?-module, called the integral clo¬ 
sure (or normalization) of R in K, and commonly denoted 
by Ok- a fc-vector space basis of K which at the same time 
generates Ok as an i?-module is called an integral basis. For 
example, when i? = Z, fc = Q, and K = Q(a) with a = 
then the canonical vector space basis {1, a, a^} of K is not 
an integral basis, because is an integral element 

of K (its minimal polynomial is — 2) but not a Z-linear 
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combination of 1, a, c?. An integral basis in this example is 
{l,a, ia^}. 

The concept of integral closure has been studied in rather 
general domains [9, 6]. To compute an integral basis for an 
algebraic number field, special algorithms have been devel¬ 
oped [7, 5]. At least two different approaches are known 
for algebraic function fields, i.e., the case when R = C[x] 
for some field C, k = C{x), and K = k[y]/{M) for some 
irreducible polynomial M € k\y\. The algorithm derived by 
Trager [10] in his thesis is an adaption of an algorithm for 
number fields, and the algorithm by van Hoeij [12] is based 
on the idea of successively canceling lower order terms of 
Puiseux series. 

The theory of algebraic functions parallels in many ways 
the theory of D-finite functions, i.e., the theory of solutions 
of linear differential operators. It is therefore natural to ask 
what corresponds to the notion of integrality in this latter 
theory. In the present paper, we propose such a definition 
and give an algorithm which computes integral bases accord¬ 
ing to this definition. Our algorithm and the arguments un¬ 
derlying its correctness are remarkably similar to van Hoeij’s 
algorithm for computing integral bases of algebraic function 
fields. 

In view of the key role that integral bases play for in¬ 
definite integration (Hermite reduction) of algebraic func¬ 
tions [10, 3, 2], we have hope that results presented below 
will help to develop new algorithms for indefinite integration 
of D-finite functions. An example pointing in this direction 
is given in the end. 

Acknowledgment. We want to thank the anonymous ref¬ 
erees for their detailed and valuable comments. 


2. INTEGRAL FUNCTIONS, 

INTEGRAL CLOSURE, AND 
INTEGRAL BASES 

Throughout this paper, let C be a computable field of 
characteristic zero, C an algebraically closed field contain¬ 
ing C (not necessarily the smallest), and x transcendental 
over C. When R is a subring of C{x), we write R[D] for the 
algebra of differential operators with coefficients in R, i.e., 
the algebra of all (formal) polynomials io-\-liD + ■ ■ ■ -\-£rD^ 
with io,... ,£r & R. This algebra is equipped with the natu¬ 
ral addition and the unique noncommutative multiplication 
respecting the commutation rules Dc = cD for all c € i? n (7 
and Dx = xD -\- 1. Typical choices of R will be C[x], C\x\, 
C{x), or C{x) in the following. 



For an operator L = Iq -\- £iD + • ■ • + irD'" £ (5[a;][Z)] 
with £r ^ 0 we denote by ord(L) = r the order of L. Recall 
that such an operator with x\ £r admits a fundamental sys¬ 
tem of formal power series solutions, i.e., the vector space 

V C (5[[a;]] consisting of all the power series / with L • / = 0 
has dimension r. When x \ £r 0, there is still a fun¬ 
damental system of generalized series solutions of the form 
exp(p(a:“^/'’))a;‘^a(a;^^®,log( 2 ;)) for some s £ N, p £ C[x], 

V £ C, a £ (7[[a;]][p]. (This notation is not meant to imply 
that a has a nonzero constant term, so the series in general 
does not start at x" but at some where i £ N is such 
that X* is the lowest order term of a.) We restrict our at¬ 
tention here to the case where p = 0 and s = 1. Moreover, 
we want to assume that v £ C (this can always be achieved 
by a suitable choice of C), to ensure that the output of our 
algorithm involves only coefficients in C. Hence we only 
consider operators L which admit a fundamental system in 
Ui/ec ®*^*5[[x]][logx]. It is well known [8] how to determine 
the first terms of a basis of such solutions for a given oper¬ 
ator L £ (7[x][II]. By a linear change of variables, the same 
techniques can also be used to find the first terms of solu¬ 
tions in Ui^gC (x — q)‘'C'[[x — a]][log(x — a)], for any given 
a £ C. More precisely, if L belongs only to C[x][-D] and 
a £ C, then these solutions are actually linear combinations 
of elements of U,/gc(® “ aYC{a)[[x — Q:]][log(x — a)]. For a 
field K with C Q K C we will use the notation 

K[[[x - a]]] := U (a: - aYK[[x - a]l[log(x - a)]. 
vec 

Observe that this is not a ring or a 7F-vector space. Also 
observe that the exponents v are restricted to the small 
field C C K, although the dependence on the choice of C 
is not reflected by the notation. We hope that the intended 
field C will always be clear from the context. 

An operator L £ C'[x][D] shall be considered integral if 
all the terms in all its series solutions remain above a cer¬ 
tain threshold. In the algebraic case, where series solutions 
involve at worst only fractional exponents, the stipulation 
of having only nonnegative exponents in all the solutions 
happens to be equivalent to the requirement that the monic 
minimal polynomial has polynomial coefficients. In the dif¬ 
ferential case, solutions involving fractional exponents cause 
factors in the leading coefficient of L regardless of whether 
the exponents are positive or negative. Therefore it doesn’t 
seem promising to use the leading coefficient of L for defin¬ 
ing integrality. Instead, we will consider the exponents of its 
series solutions. These exponents v may however belong to 
a field C which is not necessarily ordered, and there may be 
logarithmic terms. For our purposes we will use the follow¬ 
ing definition of integrality, depending on a function t which 
can be chosen according to the needs of the user. 

Definition 1. Let t:C/ZxN^'C'6ea function such that 

1. + 'L,j) £ v + 2, for every v £ C and j £ N, 

2 . l{ui +z,jl) + L(u 2 +Z,j 2 ) - + 1^2 +z, +j 2 ) > 0 

for every ui,V 2 £ C and ji,j 2 £ N, 

3. t(Z,0) = 0. 

A series f £ C[[[x — a]]] is called integral with respect to i 
if for all terms (x — a)''log(x — q)-’ occurring with a nonzero 
coefficient in f we have -\- Z,j) > 0. (For this to 

make sense the left-hand sides of the oecurring inequalities 
have to be interpreted as integers, not as elements of C.) 


The function specifies for each element u-\-Z of C/Z 
the smallest element v such that x'' log(x)-’ should be consid¬ 
ered integral. If t(iz-|-Z, j) = v, then x'' log(x)^, log(x)-’, 

... are integral and x''~^ log(x)^, x^"”^ log(x)^ ,... are not. 
The condition t(Z, 0) = 0 implies that formal Laurent series 
are integral if and only if they are in fact formal power series. 

Example 2. A natural choice for C C C is perhaps l{z -\- 
Z, 0) = a for all z £ C with 0 < K(z) < 1, and b{z-\-Z,j) = z 
for all z £ C with 0 < ^(a) < 1 when j > 1. With this 
convention, a term x’' log(x)^ is integral if and only if the 
corresponding function is bounded in a small neighborhood of 
the origin. For example, 1, x'''^, xlog(x) all are integral, 
while x~^, x'''^~^, log(x) are not. Unless otherwise stated, 
we shall always assume this choice of l in the examples given 
below. 

Proposition 3. Let a £ C and let R be the set of all C- 
linear combinations of series in (x —a)’^C'[[x —a]] [log(x —a)], 
V £ C. Then: 

1. Ln every series f £ R there are at most finitely many 
terms (x — a)^ log(x — a)^ which are not integral. 

2. The set R together with the natural addition and mul- 
tiplieation forms a ring, and {f £ R \ f is integral } 
forms a subring of R. 

Proof. 1. First consider the case when / £ (x — aYC[[x — 
Q:]][log(x — a)] for some v £ C. Let deg(/) denote the 
highest power of log(x — a) in /. Then the only possible 
non-integral terms in / can be (x — aY^'’ log(x — q)-’ for 
j £ {0,... ,deg(/)} and i £ {Q,..., b{v + Z, j) - v - 1}. 
These are finitely many. In general, if / is a linear combi¬ 
nation of some series in (x — a)‘'C'[[x — a]][log(x — a)] with 
possibly distinct v £ C, the set of all non-integral terms is 
still a finite union of finite sets of non-integral terms, and 
therefore finite. 

2. It is clear that i? is a ring. To see that the integral 
elements form a subring, let f,g £ Rhe integral. Then the 
series f -\- g cannot contain any term which is not present 
in at least one of the two summands, so all terms oi f + g 
are integral and f + g as a, whole is integral. Now consider 
multiplication: for any term (x — aY log(x — a)^ in f ■ g 
there must be some terms r in / and cr in g such that err = 
(x — Q)'^log(x — a)^, say r = (x — log(x — and 
a = (x — aY^ log(x — ay^. Since / and g are integral, we 
have /ii — t(/ii -I- Z, ji) > 0 and p ,2 — b{fj ,2 + Z,^' 2 ) > 0. The 
assumption on b in Definition 1 implies that (pi -I- P 2 ) — 
t(pi -I- p 2 + Z,ji + j 2 ) = p — t(p + Z,j) > 0. Hence all 
terms of / • p are integral, so also the product of two integral 
elements is integral. ■ 

Definition 4. Let L £ C{x)[D] and b be as in Definition 1. 

1. We call L regular if it has a fundamental system in 
C[[[x — a]]] for every a £ C. 

2. L is called (locally) integral at a with respect to b if 
it admits a fundamental system in C?[[[x — a]]] whose 
elements all are integral. 

3. L is called (globally) integral with respect to b if it is 
locally integral at a in the sense of part 2 for every 
a£C. 


Of course part 2 of this definition is independent of the 
choice of the fundamental system. In fact, L is locally inte¬ 
gral at a iff all its series solutions in x — a are integral and 
form a (7-vector space of dimension ord(I/). 

Example 5. 1. The operator {2 —x) + 2{2 — 2x + x^)D + 

4:{x — l)xD^ € Q[ 2 :][I?] is locally integral at a = 0, 
because its two linearly independent solutions 

x^ + ^x^ + ^x^ + + 0 (a;®) 

are both integral. It is also locally integral at a = 1, 
because its two linearly independent solutions 

(*-!)'/"+ 0 ((:r-l)®), 

1 - i(rr-l) + - i(rr-l)^ + 0 ((*-l)") 

are integral as well. 

The operator is also globally integral because at all a € 
C\ {0, 1} it has a fundamental system of formal power 
series, and formal power series are always integral. 

2. The operator 1 -|- xD £ Q[a;][D] is not locally integral 
af a = 0, because it has the non-integral solution i. It 
is therefore also not globally integral. 

3. The operator (—1 — 2x) -I- (x -I- 2x^)D -\- (x® -I- x'^)D^ £ 
Q[x][Z)] is not locally integral at a = 0 although all 
its series solutions are. The reason is that it has only 
one series solution in C[[[x]]] while our definition re¬ 
quires that the number of linearly independent series 
solutions must match the order of the operator. In 
other words, generalized series solutions involving ex¬ 
ponential terms, like the solution exp(i) in the present 
example, are always considered as not integral. 

Let L = £o + - ■ ■I-IrD’" £ C'[x][D] with £r ^ 0 and consider 
the quotient algebra C{x)\D]/{L), where (L) := C{x)\D]L 
denotes the left ideal generated by L in C{x)\D]. The al¬ 
gebra C{x)\D]/{L) is generated as a C'(x)-vector space by 
the basis {1, D,..., It is also a C'(x)[D]-left module, 

and we can interpret its elements as all those “functions” 
which can be reached by letting an operator P £ C{x)[D] 
act on a “generic solution” of L, very much like the elements 
of an algebraic extension field C{x)[y]/{M) can be described 
as those objects which can be reached by applying a poly¬ 
nomial P £ C{x)[y] to a “generic root” of M. A difference 
in this analogy is that in the algebraic case there are only 
finitely many roots while in the differential case we have a 
finite dimensional C-vector space of solutions. 

Definition 6 . Let L = la-\-- ■ ■-\-£rD^ £ C'[x][D] with£r yf 0 
be a regular operator and let t be as in Definition 1. 

1. An element P G A = C{x)[D]/{L) is called integral 
(with respect to l) if P ■ f is integral (with respect to i) 
for every series solution f of L. 

2. The C[x]-left module Ol of all integral elements of A 
IS called the integral closure of Clx] in A. 

3. A C{x)-vector space basis 

{Bi,...,Br}CCix)[D]/{L) 

IS called an integral basis if it also generates Ol as 
C\x\-left module. 


It is easy to see that Ol is a (^[xj-left module. Note 
however that Ol is in general not a C'[x][D]-left module, 
because the application of D may turn integral elements 
into non-integral ones (for example, D-x^^^ = when 

.(i+Z, 0 ) = i). 

Example 7. 1. The operator L = 1 — D £ Q[x][D] has 

for every a £ C one power series solution of the form 
f — 1 -b 0 (x — a). Since f is integral we have 1 £ Ol. 
Since (x — a)~^/ is not integral for any a, we have in 
fact that { 1 } is an integral basis. 

2. The operator L = 1 xD has the solution f = \. It 
is integral for every a yf 0 , but not integral at a = 0. 
However, xf = I is integral, hence x £ Ol, and in 
fact {x} is an integral basis. 

3. Whenever L has only power series solutions at every 

a G C, we clearly have {1, D,. .., C Ol . How¬ 

ever, there may still be integral elements that are not 
C[x\-linear combinations of these. For example, ob¬ 
serve that for the operator L = (x — 1) -\- D — xD^, 
which has two solutions exp(x) = l-|-x-|-|x^-|- O(x^) 
and (2x -t- 1 ) exp(—x) = x^ + 0 (x®), we have the non¬ 
trivial element i(l — D) £ Ol. Note that it is integral 
at all a ^ 0 as well. 

4 . It can also happen that 1 G Ol but D ^ Ol. For 

example, for L = (—1 -|- 2x) -b (1 — 4x)D -b 2xD^ we 
have two solutions 1 -b x -b §x^ -b O(x^) and x^^^ -b 
x 3/2 _j_ 0(x®) at a = 0. Since both are integral 

(and there are two linearly independent power series 
solutions for every a 7 ^ Oj we have 1 £ Ol. However, 
D ^ Ol, because the derivative of the second solution 
is |x“^^^-b|x^^^-b|x®^^-bO(x^), which is not integral 
since it involves the term x~^^^. An integral basis in 
this case turns out to be {l,xD}. 

5. We have produced a prototype implementation in Math- 
ematica of the algorithm described below. The code is 
available on the homepage of the first author. For the 
operator L = x^D^ -b xD — 1, it finds the integral basis 
{1, xD, xD^ — D -b i}. A fundamental system of L is 
{x,xlog(x),xlog(x)^}. 

6. Let L = 24x®D® - ISAx'^D'^ -b 373xD - 450. This op¬ 
erator has the solutions x®^^, and Our 

code finds the integral basis 

I X ’ x^ 2x® ’ X 2x2 2x® J 

In the analogy with algebraic functions, the integral op¬ 
erators from Definition 4 correspond to the monic minimal 
polynomials with coefficients in a ring, and the integral ele¬ 
ments of Definition 6 correspond to integral elements of an 
algebraic function field. Definitions 4 and 6 are obviously 
connected as follows. 

Proposition 8 . Let L £ (7[x][D] and L £ C(x)[D\ be regu¬ 
lar and assume that there exists P G C{x)[D] such that for 
every a G C we have 

{/|i-/ = 0} = {P./|L./ = 0} 

where f runs over (7[[[x — a]]] on both sides. Then P -b (L) G 
C'(x)[D]/(L) is integral in the sense of Definition 6 if and 
only if L is integral in the sense of Definition 4- 


Lemma 9. Let L = £o + • ■ • + t-rD'' G (7[x][Z)] with 7 ^ 0 
be a regular operator. Let po, ■ ■ ■ ,Pr-i G C{x) and let p = 
X — a £ C[x] be a factor of the common denominator of 
Po, • • • ,Pr-i. If Po H- \-Pr-iD^~^ G Ol then p \ tr- 

Proof. After performing a change of variables, we may as¬ 
sume that p = X. By a classical result about linear differen¬ 
tial equations (e.g., [ 8 ]), x \ Ir implies that L admits a fun¬ 
damental system bo,..., br-i in (^[[a;]] with bi = x'' + 0 (a;’^) 
for i = 0,..., r — 1. Then D^bi = i{i —1) ■ ■ ■ {i —j + l)x’'~^ + 
0(x’~~I) for i = 0,.. ., r — 1 and j = 0,..., r — 1. Let Ci 
be the largest integer such that x‘^' divides the denominator 
of Pi, let e = max{eo, •. •, fir-i}, and let i G {0,..., r — 1} be 
some index with a = e. Then piD^'bi = i\x~’^ + 
and pjD^bi = for all j 7 ^ i. Hence (po + piD + 

■ ■ ■ + pr-iD'"~^) ■ bi = i\x~’^ + ©(a;”®"*"^) is not integral be¬ 
cause —e — t(—e -I- Z, 0) = —e — t(Z, 0) = —e < 0, and hence 
po-bpi-D-l- -\-pr-iD^~^ ^ Ol- ■ 

3. ALGORITHM OUTLINE 

We shall now discuss how to construct an integral basis 
{Bo,..., Br-i} for a given regular operator L G C[a:][Z)]. 
The key observation is that van Hoeij’s algorithm for com¬ 
puting integral bases for algebraic function helds as well 
as the arguments justifying its correctness and termination 
carry over almost literally to the present setting. The re¬ 
mainder of this paper therefore follows closely the corre¬ 
sponding sections of van Hoeij’s paper. 

The algorithm computes the basis elements Bo,..., B ,— 1 
in order, at each stage d G {0,..., r — 1} starting with an ini¬ 
tial conservative guess for Bd and refining it repeatedly until 
an operator Bd is found which together with Bo,..., Bd-i 
generates the C[x\-\eit module consisting of all the elements 
of Ol corresponding to operators of order d or less. Al¬ 
though parts of the calculation take place in the large field C, 
it will be shown that the elements Bi in the resulting inte¬ 
gral basis always have coefficients in the small field C, in 
which the coefficients of the input operator L live. 

It is not hard to find a suitable Bq: For each root a £ C 
of the leading coefficient £r of L, compute the first terms of 
a basis {for,... , 6 r} of solutions in C[[[a: — a]]]. Determine 
the smallest integer Ca such that {x — aY°‘bi is integral for 
every i according to the chosen l. Then Bo can be set to the 
product of {x — a)®" over all a. Since Ca = ea whenever a 
is a conjugate of a, it follows that Bo belongs to C{x). 

The outline of the algorithm is now given on a conceptual 
level. In Section 5 a more detailed description of steps 5-7 
will be given. 

Algorithm 10. 

INPUT: A regular operator L = lo + irD^ £ C'[ 2 ;][D] 
with £r ^ 0 

OUTPUT: {Bo,..., Br-i} C C{x)[D\/{L), an integral basis 
of C{x)[D]/{L). 

1 Set s to the squarefree part of £r. 

2 Set Bq to the zero-order operator described above. 

3 For d = 1,..., r — 1, do the following: 

4 Set Bd = sD Bd-i. (Also Bd = s'^D'^Bo would work.) 
Consider 

E = {A£Ol :oTd{A) <d}\{C[x]Bo + --- + C[x]Bd). 

5 While B 7 ^ 0, do the following: 


6 Construct A G B of the form 

A = — (floBo -|- • • • -|- Od-iBd-i + Bd) 

with ao,. .., ad-i,p G C[x\. 

7 We have 

C[x]Bo -b • • • -b C[x]Bd-i + C[x]Bd 
C C{x\Bo -b • • • -b C[x\Bd-i -b C[x\A C Ol- 

Replace Bd by A. (This makes B strictly smaller.) 

8 Return (Bo,..., Br-i}. 

In the rehned version of the algorithm, we will see that 
the set B is never explicitly constructed. Instead, it suffices 
to be able to solve the following subproblems. First, we need 
to decide whether B = 0 for recognizing the termination of 
the loop in lines 5-7; this is discussed in Section 5. Second, 
we need to show the existence of an element A G B of the 
form required in step 6 whenever B 7 b 0; see Section 4. In 
Section 5 we explain how such an A is constructed. Finally, 
the termination of the loop in lines 5-7 is proved in Section 6 . 
Except for these issues, the correctness of the algorithm is 
obvious. 

4. EXISTENCE OF A IF / 0 

The arguments in this section are almost identical to those 
in [12]. Nevertheless, for sake of completeness, we formulate 
them here for the differential case. 

In the d-th iteration of the algorithm we can assume by 
induction that Bq, ..., Bd-i with ord(Bi) = i for all i form a 
(7[a;]-left module basis of all integral elements of order up to 
d — 1 . We consider the case where the current choice of Bd, 
together with Bq, ..., B^-i, does not generate all integral 
elements of order up to d, i.e., E ^ <!). Recall that 

B = { A G Oi : ord(A) <d}\ {C\x\Bo -b ■ • ■ -b C\x\Bd). 

We need to show that there exists an integral element A £ E 
which can be written in the form i (aoBo -b • • • -b aaBd) with 
ao,... ,ad,p G C[x] and = 1. The idea is as follows: 
starting from an arbitrary element A G B, we construct, in 
several steps, simpler elements in B until we obtain one with 
the desired properties. 

Lemma 11. //B 7 b 0, then there exists A £ E of the form 
A = ——— (aoBo -b • • ■ -b Od-iBd-i -b OdBd) (1) 
with a £ C, ao,..., ad-i,ad G C[x\. 

Proof. Let A £ E, say A = aoBo -b • ■ • -b adBd for some 
Oi £ C{x). Since A ^ C[x\Bo -b • • • -b C[x\Bd, at least one 
Oi must be in C{x) \ Clx]. Let p G C[x\ be the common 
denominator of all the ai, and let a G C be a root of p. 
Then ^^A has the required form. To see that it belongs 
to B, notice that G C[x\ and Ol is a ( 7 [ 2 :]-module, and 
that ^^A ^ Cla^lBo H--b C[x\Bd. ■ 

Lemma 12. If A £ E and P G C[x\Bo -b ■ • • -b C[x\Bd, then 
A-bP G B. 

Proof. A £ E £- Ol and P £ C'[a;]Bo -b • ■ ■ -b C[x]Bd C Ol 
implies that A-bP G Ol. It is also clear that ord(A-bP) < d, 
because ord(A) < d and ord(P) < d. Finally, to show that 



A+P ^ C[x\Bq + - ■ ■ + C[x\Bd, assume otherwise. Then also 
A = {A-\- P) — P £ C[x\Bq + ■ ■ • + C[x\Bd in contradiction 
to A £ E. ■ 

Lemma 13. If E contains an element of the form (1), then 
it also contains such an element with ao,.. ■ ,ad-i £ C and 
ad = 1 . 

Proof. Let A = (aoBo + • ■ • + OdBd) G L be of the 
form (1). For each i = 0,..., d, write ai = {x — a)pi + a[ 
with Pi £ C\x\ and a' £ C. By Lemma 12, A £ E implies 
A £ E tor A := -^A— ^aoBo-£ ■ ■ ■ Pad—iBd—i PadEd) ■ Since 
Bq, ..., Bd-i are assumed to generate the submodule of all 
the elements of Ol of order at most d — 1, we have a'd 7^ 0. 
Dividing A' by a'd yields an element of E of the requested 
form. ■ 


Lemma 14. If E contains an element of the form (1) with 
ao, ■ ■ ■ ,ad-i £ C and = 1, then it also contains such an 
element with ao,, Od-i £ C{a) and = 1. 

Proof. Let A £ E he ot the form (1) with ao,..., ttd-i £ C 
and = 1. Since C is necessarily a C'(Q;)-vector space, there 
are some C'(a)-linearly independent elements eo,...,e„ of 
C such that ao,... ,ad all belong to F = eoC'(a) + • ■ • + 
e„C(a). We may assume eo = 1. Consider a fundamental 
system b\,... ,br £ (7(0?) [[[x — a]]] of L. Then each A ■ bj 
has coefficients in V and, since A £ E £ Ol, only involves 
integral terms. For an element n € F let us write [ei]v 
for the coordinate of v with respect to Ci. By the linear 
independence of the Ci over C{a), the series [ei](A ■ bj'j = 
{[ei]A)-bj obtained from A-bj by replacing each coefficient by 
its Ci-coordinate will be integral. In particular, the operator 
Ao = [eo]A £ C(a)[x][D] must belong to E. Because of 
[eo\ad = [eo]l = 1, it meets all the requirements. ■ 

Lemma 15. If E contains an element of the form (1) with 
ao,... ,ad-i £ C{a) and ad = 1, then it also contains such 
an element with ao,..., ad-i £ C[x\ and = 1. 

Proof. For every n > 0 we have x — a \ x’' — a’' in C[x], 
and thus also x — a \ p{x) — p{a) for p £ C[x] \ C. There¬ 
fore, if we view the at £ C{a) as polynomials in a, then 
replacing a in them by x amounts to adding some polyno¬ 
mial multiple of {x — a) to them. This change means for 
A = -i^^iaoBo -I- • ■ • -I- ad-iBd-i + Bd) that adding a suit¬ 
able element P £ C{q)[x]Bo + - ■ ■ + C{a)[x]Bd-i £ Ol turns 
A into an operator of the requested form. By Lemma 12, 
this new operator also belongs to E. m 


Theorem 16. If E ^ then there exists an element A £ E 
of the form 


A = — (aoBo -I- ■ • ■ -I- ad-iBd-i 

p 


+ Bd) 


with p £ C[x\ an irreducible factor of Ir and ao,... ,ad-i £ 
C[x] such that deg(ai) < deg(p) for all i. 


Proof. The assumption i? 7^ 0 in combination with Lem¬ 
mas 11, 13, 14, and 15 implies that E contains an element of 
the form (1) with ao,..., ad-i £ C\x\ and ad = 1. Further¬ 
more, Lemma 9 implies that a is a root of ir. Let p\ Ir he 
the minimal polynomial of a. We claim that A ^B £ E 
where B ;= aoBo H--I- ad-iBd-i + Bd. 


To prove this, we have to show that for every a £ C 
and every solution b £ (7(0;) [[[x — o]]] of L we still have 
that ^ ■ fo is integral. When d is not a root of p, this is 
clear because 1/p admits an expansion in C[[x — d]], and 
multiplication of the integral series B ■ b hy a. formal power 
series preserves integrality by Proposition 3. When d = o, 
write p = (x — a)q for some q £ (7[x] with x — o f q and note 
that 1/q admits an expansion in C[[x — o]] and y;^B ■ b is 
integral, so ^B ■ b is integral too. When d is a conjugate 
of a, note that • b must be integral, because if it were 

not, then for the series b £ (7(q)[[[x — o]]] obtained from b 
via the conjugation map that sends d to q we would have 
that ■ b is also not integral, in contradiction to our 

choice of ao,..., ad. Therefore the same argument as in the 
case a = a applies. 

This completes the proof of the claim. To complete the 
proof of the theorem, note that the claimed degree bounds 
on Oi can be ensured by Lemma 12. ■ 

5. CONSTRUCTION OF A IN STEP 6 

In the previous section we have demonstrated that in 
step 6 of the algorithm it suffices to search for an integral 
element A of the form 

^ ~ [aoBo -l- ■ • ■ ad—iBd—i + Bd) 

where ao,... ,ad-i,p £ C[x], p \ ir and deg(ai) < deg(p). 
Conversely, this means that if no such A exists, the set E is 
empty. 

For each irreducible factor p of ir one can set up an ansatz 
for A with undetermined coefficients ao, ■.., a^-i. We want 
to find ao,..., Od-i such that A ■ f is integral for all solu¬ 
tions / of L. Note that we need to enforce integrality only 
for series solutions in x — a where a is a root of p. Choosing 
a fundamental system hi,... ,hr of such solutions, comput¬ 
ing the first terms of Bj ■ bi, plugging them into the ansatz, 
and equating the coefficients of all non-integral terms to zero 
yields a linear system for ao,... ,ad-i. If this system does 
not admit a solution, one knows that no such A with de¬ 
nominator p exists. 

In summary, the loop in lines 5-7 of Algorithm 10 can be 
described in more detail as follows. 

5a Let Q £ C he a set containing exactly one root a £ C 
for each irreducible factor p of ir. 

5b While (3 7^ 0, do the following: 

5c For all a G (3, do the following: 

6a Let bi,... ,br be a fundamental system of L in 

C(a)[[[x-«]]]. 

6b With variables ao,..., Od-i, form the series 

(aoBo • -l- ad—iBd—i A Bd)hi 
for i = 1 ,..., r. 

6c Construct a linear system for ao,..., ad-i by equat¬ 

ing the coefficients of all the non-integral terms in 
these series to zero. 

7a If the system has a solution (ao,..., Od-i) £ C(a)'^: 

7b Let p be the minimal polynomial of a over (7. 

7c Replace each ai £ C(a) = (7[x]/(p) by the corre¬ 

sponding polynomial in (7[x] of degree less than 
deg(p). 


7d Replace Bd by i(ao-Bo H-h ad-iBd-i + Bo). 

7e Otherwise 

7f discard a from Q. 

Despite being more detailed than the listing given in Al¬ 
gorithm 10, these lines are still somewhat conceptual. An 
actual implementation cannot just “let” bi be some infinite 
series object, and it does not need to. What we need are only 
the terms of bi that give rise to some non-integral terms of 
(aof3o + ' • ■+ad-iBd-i+Bd)bi. These are only finitely many 
by Proposition 3. In Section 7 we address the question how 
many terms of bi we need to compute. 


6. TERMINATION 

The termination of van Hoeij’s algorithm [12] is estab¬ 
lished by the observation that the degree of a certain poly¬ 
nomial, starting with the discriminant ReSy (M, ^■), de¬ 
creases in each iteration of the main loop. In the case of 
D-finite functions, the role of the discriminant is played by 
the Wronskian and a generalized version of it. Recall that 
the Wronskian of the functions /i (x),..., fr{x) is dehned as 
the determinant 


fl{x) 

h{x) 

fr{x) 

fl{x) 

Mx) 

frix) 

t-^\x) 


■■ 


Definition 17. Let L G C'[x][D] be regular and letbi,... ,br 
be a fundamental system of L in (7][[x — a]]] for some a £ C. 
For Bq,...,B, — 1 G C{x)[D]/{L) we define the generalized 
Wronskian at a, as 


VfrL,a{Bo, . . . , Br-l) 


Bo ■ bi 


Br — l • 


Bq ' br 



Note that the generalized Wronskian 'NiL,a{Bo, ■ ■ ■, Br-i) 
belongs to (7[[[x — a]]] and that the choice of a different 
fundamental system instead of bi,...,br only changes its 
value by a nonzero multiplicative constant, which will be 
irrelevant for our purpose. 

For the special choice Bi — D', the generalized Wron¬ 
skian wTL,a{i, D,... , D' reduces to the Wronskian (2) 
with fi = bi. It is well-known and easy to check that the 
classical Wronskian (2) of 6i,... , br satisfies the hrst-order 
equation IrD W +£,— iW = 0 and hence is hyperexponential. 
Since the generalized Wronskian can be obtained from the 
usual Wronskian by elementary row operations over C'(x), 
it is clear that also the generalized Wronskian is hyperexpo¬ 
nential. 


Theorem 18. Algorithm 10 terminates. 


Proof. First observe that during the whole execution of the 
algorithm, Bo,...,Br-i G C{x)[D]/{L) are integral, i.e., 
Bq ■ f,...,Br-i ■ f are integral for any series solution / 
of L according to Definition 6. (Actually, the Bd’s are 
constructed one after the other, but they can be initialized 
with Bd = s’^D'^’Bq.) This means that, at any time and for 
any a £ C, the generalized Wronskian wrL,Q,(Bo,.. ■, Br-i) 
is integral, as it is the sum of products of integral series 
(see Proposition 3). Since it is hyperexponential, it follows 
that it has no logarithmic terms. Every nonzero term of 


wri:,,ct(Bo, ■ • ■, B ,— i) is therefore of the form (x — q)'' with 
fi = t(/r -I- Z, 0) + m for some nonnegative integer m. For 
each a G (7 let ma be the smallest such integer. Now let n = 
^asQ where Q is defined as in step 5a. Each time Bd 
is updated in the algorithm (either in step 4 or in step 7d), 
none of the ma can increase and exactly one of them strictly 
decreases, so also n decreases. More precisely, if for exam¬ 
ple Bd is replaced by ^{aoBo -I- • ■ • -I- Od-iBd-i + Bd) in 
step 7, then wTr.aiBo ,..., Bd) is divided by p (recall that 
p is a non-constant polynomial in C[x\). But the ma can¬ 
not become negative as this would violate the integrality of 
'^^L.ai.Bo,... ,B ,— i). Therefore the algorithm must termi¬ 
nate. ■ 


7. BOUNDS 

In the algebraic case, van Hoeij [12] gave a-priori bounds 
on the orders to which the bi have to be calculated. His 
algorithm computes their terms at the very beginning once 
and for all in order to avoid their recomputation inside the 
loop. He also suggested that the terms of Bj ■ bi for j < d 
should not be recomputed but cached. 

Nowadays, in an object-oriented programming environ¬ 
ment, the algorithm can be implemented in such a way that 
recomputations of series terms are avoided even when no 
a-priori bound on the truncation order is available, via the 
paradigm of lazy series [4, 11]. 

Nevertheless it is desirable to have a-priori bounds avail¬ 
able also in the D-finite case. A rough bound follows imme¬ 
diately from the discussion in Section 6: as we have seen, 
the Wronskian wTr.a (Bq, sDBq, . .., s^~^B q) gives a 
denominator bound for the elements of the integral basis. 
More refined bounds are elaborated in the following. 

Let a G (7 be a root of the leading coefficient and 
{bi ,... , br} C C'(a)]]]x — a]]] be a fundamental system of L: 

OO 

6i = ^ 6i,fe(log(x - a)) (x - a)"*"'"*, 6i.o / 0, (3) 

k = 0 

where bi^k € C(a)]log(x — a)] are polynomials in log(x — a) 
such that for each i the degrees of bi.o, 5i,i,.. . are bounded 
by some integer di. According to step 5c, we have to consider 
each a £ Q separately, so for the rest of this section we fix 
such an a. 

In step 6a we want to replace bi,... ,brhy truncated series 
ti,... ,tr of the form 

ti = ^ fci.fe (log(x — a)) (x — a)''''^^ with Ni £ N. (4) 

k=0 

The bounds Ni must be chosen such that this replacement 
does not change the result of the algorithm. The only criti¬ 
cal step is when bi,... ,br are used to test the integrality of 
certain elements from the algebra C{x)[D]/{L), which are 
not known in advance. Theorem 20 below gives a sufficient 
condition that allows us to use ti instead of bi in the inte¬ 
grality test, by asserting that its answer does not change, 
whatever element of C{x)[D]/{L) we consider. For brevity, 
let R denote the ring C(q)]]x — a]]]log(x — a)] in the subse¬ 
quent reasoning. 

Lemma 19. Let {6i,...,br} C C(q)]]]x — q]]] be a fun¬ 
damental system of the form (3) with Vi as above, and let 
Wb = {D^ ■ bi)i<i<r,o<j<r- Then we can find an m G N sueh 






that 


det(W"6) = Y,^k{x- a)-i+-+---(’-i)/2+-+'= 

k = 0 

With Wo ^ 0. 

Proof. For the (i, ji)-entry of Wb we have 

{Wb)i,i = -bi^ix- ar-^+^R 

and therefore 

det(VF6) £ (x- 

Note that det(Wj,) 0 because it is precisely the Wronskian 
of &i,..., 6r. It follows that a unique m > 0 with the desired 
property exists. ■ 

Theorem 20. Let L £ C{x)[D] be an operator of order r 
and {bi,..., br} C C'(a)[[[x — a]]] be a fundamental system 
of L with Vi and di as above. Moreover, let m £ N be as in 
Lemma 19 and let Ni,... , Nr G N be given by 

Ni = m + max ( b{vi — Vj + Z, fc) — {vi — Vj) ). 

l<j<r V / 

0<fc<dj+r 


of W^, ^ is r — 1. On the other hand, it is easy to see that 
Zi^j £ (a: — R, so it follows that 

. £ (rr - aY^-'^^+^'-^+^R, (8) 

and that herein log( 2 ; — a) appears with exponent at most 
di + r — 1. By our choice of Ni the series in (8) is integral 
for all 1 < i, j < r and therefore the whole matrix . 

Part 2: Now assume that B - b is not integral. Then from 

B-t^ (Id, - ZWY^) {B-b) = B-b - {ZWY^) {B ■ b) 

it follows that B ■ t is non-integral as well. To see this, 
let n be the largest integer such that a term of the form 
(x — Q/)dM+z,fe)-n — ctY appears in B • 6 for some /r £ C 
and fc £ N. Let i be an index such that a term of the given 
form appears in B • 6i with nonzero coefficient. This term 
cannot be canceled in 

r 

B-ti = B-bi-J2 {ZWb-^)^^^{B-b^) 

because all terms of the series (ZW^^)^ ^ are of the form 

(x — log(x — aY with t! > 1 by our choice 

of Ni. So also B ■ t is not integral. ■ 


If ti is the truncation (4) of bi at order Ni, for 1 < i < r, 
then for all B £ C{x)[D]/(L) we have the equivalence: 

yi: B ■ bi is integral <==> yi: B ■ ti is integral. (5) 


Proof. We introduce the matrix 144 = {D^ ■ fei)i<i< 7 -,o<j<r 
as before, and the short notation B ■ b — {B ■ bi,..., B ■ br). 
Analogously we define 144 and B ■ t. A vector resp. matrix 
is called integral if all its entries are integral. If c is the 

coefficient vector of B, i.e., c • (1, D,..., D’ = B, then 

we have B - b — 144 c and B -t = 144 c. Combining these two 
equations we get 

B ■t = WtWY^ {B-b). ( 6 ) 

Setting Z = 144 — 144 yields 

WtWY^ =ldr- ZWY\ ( 7 ) 


The proof is split into two parts, according to the two 
directions of the equivalence (5). 

Part 1: If we assume that B-b is integral, then (6) exhibits 
that the integrality of WtWY^ is a sufficient condition to 
conclude that also B ■ t is integral, using Proposition 3. By 
(7) it suffices to show that ZW^^ is integral. First of all we 
have to argue that £ (^((^^[[[x — since otherwise 

Definition 1 would not be applicable. In Section 6 we have 
remarked that the Wronskian det(114) is hyperexponential. 
In particular, it involves no logarithmic terms and therefore 
is invertible in (7(0;)[[[x — o]]]. Using Cramer’s rule we find 
that 


^ detWb 


R, 


where is the matrix obtained by deleting row j and 

column i from Wb. So the entries of are series in 

(7(o)[[[x — o]]]. The fact that detlFjJ^’*' satisfies a differ¬ 
ential equation of order less than or equal to r implies that 
the highest power of log(x —o) that can appear in the entries 


8. COMPARISON WITH THE 
ALGEBRAIC CASE 

We have shown that the underlying ideas of van Hoeij’s 
algorithm for computing integral bases of algebraic function 
fields apply in a more general context. Indeed, it is fair to 
regard van Hoeij’s algorithm as a special case of our algo¬ 
rithm, since every algebraic function is also D-finite. Recall 
that an algebraic function field C{x)[y]/(M) with some ir¬ 
reducible polynomial M of degree d becomes a differential 
field if we set D ■ c = 0 for all c £ (7, D • x = 1, and 

Y-M 

D ■ y := — mod M. 

Y-M 

dy 

Since C{x)[y]/{M) is also a (7(x)-vector space of dimen¬ 
sion d, it is clear that any d -|- 1 elements must be C(x)- 
linearly dependent. This implies the existence of an opera¬ 
tor L £ C{x)[D\ of order at most d with L ■ y — 0. Usually 
there is no such operator of lower order, which means that 
y,D-y,..., D‘^~^ ■ y are (7(x)-linearly independent and thus 
a basis of C{x)[y]/{M). In this case, a vector space basis 
{Bi,..., Bd} U C(x)[t/]/ (L) is an integral basis in the sense 
of Definition 6 if and only if {Bi-y,..., Bd-y} Y C(x)[j/]/(M) 
is an integral basis of the algebraic function field in the clas¬ 
sical sense. 

When y £ C(x)[y\/{M) is annihilated by an operator L 
of order less than d, we can compute the minimal-order op¬ 
erators Lo,..., Ld-i which annihilate y^,..., y‘^~^, respec¬ 
tively, and take L = lclm{Lo,..., Ld-i). Then the (7(x)- 
vector space generated by all solutions of L is the whole 
field C{x)[y]/{M), and if {Bi,..., Bn} is an integral basis 
for L, then {Bi -y^ : i = 1,..., n, j = Q,... ,d—l} generates 
the (7[x]-module of all integral elements of (7(x)[j/]/(M). 

As a less brutal approach, we can simply replace y by 
some other generator of the field. In practice, most field 
generators will have an annihilating operator of order d, but 
none of smaller order. 




Example 21. An integral basis for the field Q{x)[y]/{M) 
with M = is {l,y, The lowest-order differ¬ 

ential operator annihilating y is L = 3xD — 2, which is not 
useful because its order is less than the degree of M. 

Instead, let us try Z = 1 + y y^ as generator. We have 
Q{x)[y]/{M) = Q{x)[Z]/{N), where N = Z^-3Z‘^- 3(x^ - 
1)Z — a;'* + 2a;^ — 1 is the minimal polynomial of Z. Given 
N instead of M as input, van Hoeij’s algorithm finds the 
following integral basis for Q(a:)[Z]/(A''); 

^ ^ _ + 2)Z 1 

' x{x — l){x1) x{x—l){x + l) X 

The lowest order annihilating operator of Z is L = 9x^D^ + 
9xD^ — D. It has the right order and our Mathematica im¬ 
plementation returns the integral basis {l,xD,xD^ + |-D}. 
We can express its derivatives as polynomials in Z, using 

^ + 2(2x^ + l)Z 

3x{x — l)(x + 1) 

and obtain the following integral basis for Q(x)[Z]/(Ai); 

-2Z^ + 2(2a:^ + 1)Z 8{-Z^ + {x^ + 2)Z + - 1) ) 

3(a; — l)(a; + 1) ’ 9a:(x — l)(a; + 1) / ’ 

Applying a change of basis with the unimodular matrix 




1 

8 


-12 

0 

0 



gives the integral basis (9) computed by Maple. 


One of the features of integral bases for algebraic function 
fields is that they allow an extension of the classical Hermite 
reduction for integration of rational functions to the case 
of algebraic functions. This was observed by Trager [10]. 
In order to make this work, Trager requires that both the 
integral basis as well as the integrand should be “normal at 
infinity”. This corresponds to the condition in the rational 
case that the rational function to be integrated must not 
have a polynomial part. Trager shows that normality of 
the integrand can always be achieved by applying a suitable 
change of variables, and he gives an algorithm that turns an 
arbitrary integral basis into one that is normal at infinity. 
After that, the Hermite reduction process looks very similar 
to the rational case. We give here an example for a non- 
algebraic D-finite function. 

Example 22. Let L = (22; +1) — (42;^ + l)D-\-2{2x—l)xD^; 
its solutions are exp( 2 :) and y/x, but we will not use this in¬ 
formation. Let us just write y for a solution of L. An inte¬ 
gral basis ofOr is given by {1, ^^^{2xD — 1)}. Let ujq := y 
and cui ;= 2 x^-i i^xD — 1) ■ y and consider the function 

_ aoujo + aicui 

J ^ 

where ao = 42;^ + 37x — 11, ai = —28x^ + 402:^ — x — 1, 
u = A, V = (x — l)x, m = 2. 

Hermite reduction consists in finding 6o,6i,co,Ci G Q[ 2 :] 
with 

aoujo + aiuji _ / bowo + biuui N ' ^ cowo + ciuui 
UV'^ V y-m-l J yym-l 

After working out the differentiation, multiplying by ww"*, 
and taking the whole equation mod v we are left with the 


constraint 

aouJo + aicJi = bouv”^ ( ) + biuv^ ( ) mod u 

\ym-lj \ym-l ) 

For the derivatives of wo and wi we have 

^ 1 l-2x 

Dluo = —Wo -T- wi, Dwi = wi, 

2x 2x 

so that the previous constraint can be rewritten to 

aowo + aiwi = —^bou{3wo + wi) — 2biuwi mod v. 

Plugging in oo,ai and u and comparing coefficients of Wi 
leads to the linear system 

fAlx - 11\ {2-&X 2 - 22;\ {bo\ 

[llx-l) = [ 0 

which has the solution bo — ^(42; + ll), for = |(22: —1). Next 
we find that 

/ bowo + biwi \ ' _ cqWo + CiWi 

J \ ym-l ) ~ uym-l 

for Co = 0, Cl — 0. Consequently, we have found that 

r _ (ll + 42:)aJo + 5(22:-l)a)i _ 5 , _ 22: + 3 

J 8{l—xyx'^ X — 1^ {x — l)x^ 

The same answer could have been found using an algorithm 
of Abramov and van Hoeij [1], using a different approach. 
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